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Abstract 

In this study, it is theoretically proven that the expected value of maximum loss of frac- 
tional Brownian motion (fBm) up to time 1 with Hurst parameter [1/2, 1) is bounded 
above by 2/y/rr and below by 1/ ' \[tt. This result is generalized for fBm with H £ 
[1 /2, 1) up to any fixed time, t. This also leads us to the bounds related to the distribu- 
tion of maximum loss of fBm. As numerical study some lower bounds on the expected 
value of maximum loss of fBm up to time 1 are obtained by discretization. Simula- 
tion study is conducted with Cholesky method. Finally, comparison of the established 
bounds with simulation results is given. 

Keywords: Cholesky Method, Hurst Parameter, Fractinal Brownian Motion, 
Maximum Loss, Sudakov-Fernique Inequality. 



1. Introduction. 

Let H be a constant in the interval (0, 1). Fractional Brownian motion (fBm), 
-Bj , t > with Hurst parameter H, is a continuous, centered Gaussian process with 
covariance function 

E[B?B?] = l(t" I + t ? B -\t-s\™) (1) 

For H — 1/2, fBm corresponds to standard Brownian motion. 
A standard fBm, B H , t > has the following properties, JH: 

• = and E[B?\ = for all t > 0. 

• B H has homogenous increments, that is B^ +s — has the same law as Bf , 
fors,f > 0. 

• B H is a Gaussian process and E[(B? ) 2 ] = t 2H , t > 0, for all H £ (0, 1). 

• B H has continuous trajectories. 
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1.1. Correlation between two increments and long range dependence 

For H = 1/2, the process Bf, t > corresponds to a standard Brownian mo- 
tion, in which the increments are independent. By the definition of fBm, we know the 
covariance between B H (t + h) — B H (t) with s + h <t and t — s = nh is 

PH {n) = h 2H [(n + l) 2H + (n - 1) 2H 2n 2H ], (2) 

We observe that two increments of the form B H (t + h) - B H (t) and B H (t + 2h) - 
B H (t + h) are positively correlated for H > 1/2, and they are negatively correlated 
fori? < 1/2, Jl. 



Definition 1. A stationary squence (X n ) ne m exhibits long-range dependence if the 
autocovariance function p(n) := cov(Xk, Xk+ n ) satisfy 



lim 



P( n ) _ 1 



for some constant c and a £ (0, 1). In this case the dependence between Xk and Xk+ n 
decays slowly as n tends to infinity and 



n=l 

Hence, it can be obtained that the increment Xk — Bf (k) — Bf (k — 1) of Bf 
and Xk+n — Bf (k + n) — B^(k + n — 1) of Bf have the long-range dependence 
property for H > 1/2 since 

p H (n) = I[( n + lfH + (n- i)™ - 2n 2H ] » H(2H - l)n 2H - 2 

as n goes to infinity, yfl. In particular 

lim ; r orr ^ = 1. 

n^ooH(2H - l)n 2H - 2 

1.2, Self similarity property 

FBm possesses the self-similarity property, that is for any constant c > 0, 

(B? t ) t > l ^(c H B?) t > . (3) 

1.3. The model 

Fractional Brownian motion (fBm) is used in modelling many situations such as 
modelling temperature, characters of solar activity, values of the log returns of a stock, 
price of electricity, fl. FBm with H € [4, 1) is used in finance due to its long-range 
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dependence property, fl. The Black-Scholes model for the values of returns of an 
asset using fBm is given as, 

Y t = Y a exp((r + fi)t + aB?) < t < T, (4) 

where Yq is initial value, r constant interest rate, /j, constant drift and a constant dif- 
fusion coefficient of fBm. Here, fBm is denoted by Bf, t > 0. Black-Scholes model 
can also be constructed with Markov processes such as Brownian motion. The ad- 
vantage of modeling with fBm to Markov proccesses is its capability of exposing the 
dependence between returns. Real life data for a volatile asset, display long-range de- 
pendence property. For this reason, modelling with fBm is more realistic compared to 
Markov processes. 

In financial markets, investors would be interested in the maximum possible loss as 
a measure of their highest possible risk. Some bounds on the expected value of maxi- 
mum loss of fBm already exists, [6]. There are only bounds for fBm since exact results 
are not possible to find. Our aim is to improve these bounds. 

In this study, we provide a new and closer theoretical upper and lower bounds on 
the expected value of maximum loss. We obtain some numerical lower bounds on 
the expected value of maximum loss using discretization and order statistics of multi- 
variate normal variables. Simulation results are obtained by Cholesky method. As a 
conclusion, we provide comparison of the established bounds. 

1.4. Notation and Preliminaries 

Let Bf, t > be a fBm defined on the probability space (il, J 7 , P). Let 

Sf := sup B^ be the supremum, if :~ inf B^ be the infimum of fBm, — 

o<v<t * a<v<t 

Sf — if be the range, the difference between the supremum and the infimum, Mf := 
sup (S„ — B^) be the maximum loss of fBm. Some theoretical bounds found on 

0<u<v<t 

the expected value of supremum and the expected value of maximum loss of fBm, 

• An upper bound on the expected value of supremum was given as, 11811 : 

E(Sf) < ^Lt H . (5) 

• Additionally, a lower bound on the expected value of infimum and an upper 
bound on the expected value of range are,[0]: 



-^L? < E (I t H ), E(R?) < ^t H . 



(6) 



In II12I1 using Sudakov-Fernique inequality upper and lower bounds on the ex- 



pected value of supremum of fBm were given as 
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• Combining all the above and the fact that E(S^) equals —E(I^) provided 
bounds on the expected value of maximum loss of fBm up to time t, Jg]: 

^L* < E(M t H ) < ^t H . (8) 



2. Main Results 

In this section, we work on improving the bounds given in equation ([S]). We prove 
that the expected value of maximum loss of fBm up to time t is bounded above by 
2t if /y / 7r which is a closer bound than 2\p2t H j ^JH. Also, we prove that the expected 
value of maximum loss of fBm up to time t is theoretically bounded below by t H /^/tt 
which is again a closer bound than \/2t H /2^/tt. Our second new result is to provide 
some numerical lower bounds on the expected value of maximum loss of fBm. Finally, 
we give the results of our simulation study in order to be able to compare the established 
bounds. 



2.1. Theoretical Results 

Here, we prove that the expected value of maximum loss of fBm up to time 1 with 
He [1/2, 1) is bounded above by 2/y/ir and it is bounded below by l/y/n. These 
new bounds provide closer and improved values than the bounds given in equation ([8). 

Theorem 1. For fBm up to time 1, with Hurst parameter 1/2 < H < 1, 
-==E{ sup V2Bl) < E(M?) < E( sup V2B] /2 ) = 

V 71 " 0<s<l 0<s<l V 7T 

Proof. Now, {B^ — B^} are centered Gaussian processes on seperable space [0, 1] 
and 

E{{B*-B»)-{B*-B»)f = E((B»-B»f) 

2E((B*-BX)(BZ-B*)) 

+ E{{B»-B*f) 

i \ 2H i i ' '\ 2H i i n 2H 
= \u — v\ + \u — v | + \u — u I 

— \u — i/| 2ff — \v — u'\ 2H + \v — v'\ 2H (9) 



For all u, v G [0, 1] and 1/2 < H < 1 note that, 



\u — v\ 2H < \u — v\ = u + v — 2min(u, v). (10) 

Considering all 6 cases such as < u < v < u' < v' < 1, < u < u' < v < v' < 1, 

0<«<«'<t)'<»<l, 0<u'<v'<u<v<l, 0< u' <u<v'<v<l, 
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< u' < u < v < v' < I , one by one we see that, 



eub%-b*)-{b5-b5)Y 



< 



< 



u — v\ + \u — v I + \u — u | 
U — 7J | — |u — u I + |u — v I 



,2H 



|2H . 
14 — V\ + | 



m - v \ 



n2H I /|2// 

u-ti | +\v — v\ 



\u — v\ + 1 14 — t> | 
+ |l4- 1i'| + |l> - l/| 

< 2 max(ii, u, u , v ) + 2 min(u, u, u , i> ) 
— 4mm(min(w, u, it , u ), max(it, u, u , 1! )) 

(11) 

And, for all u, w € [0,1], note that 

E{^2B l J 2 - V2By 2 ) 2 = E(2(B 1 J 2 ) 2 ) + E{2(B l J 2 ) 2 ) - 2E(2B 1 J 2 B 1 / 2 ) 

= 2 | it — i>| = 2u + 2« — 4min(u, v) (12) 

Therefore, we have 

E{{Bu — By) — {B^i — £v )) 2 < 2max(it,i;,i4 / ,i; / )+2min(i4,i;,it',t; / ) 

— 4mm(min(u, v, u , v ), max(u, v, it', v )) 
= EiV^B 1 

V v I] 

Hence, by the Sudakov-Fernique inequality given in 



max{it.y,u'.w'} 



V2i? 1/2 , , ,,) 2 



E(M?) < E{ sup V2B 1 / 2 ) = E(V2Sl /2 ) = -== 

0<s<l V 71 " 



(13) 



is satisfied. 

Next, we prove lower bound of the expected value of the maximum loss of ffim up 
time 1 is 1/ y/n. From equation ©, we have 

n 2H 



E((BZ~B*)-(B%-B%)) 2 = \u-v\ 



2H 



I / / \2H 
It — VI 



I f\2H I f\2H I f\2H 

— 1 14 — If | — \V — U | + |t! — U | 



Suppose 0<i4<u<i4'<i/<lis satisfied. We would like to show that 

E{{BS - flf ) - (2# - i^)) 2 > [14 - if" + |i/ - v'\ 2H 
Equivalently, our aim is to show 



Now, 



|14 


-u> 


2H + 


|«- 


i\2H 

-v\ - 


14- 


-v\ - 


\v-u'\ 2H > 0. 


[14 


-u 




|« 


n 2H 

-v\ - 


[14 


n 2H 

-v\ - 


|« - u'\ 2H 




-u 




| V 


-v\ - 


[14 


-u' + u' 


-v'\ 2H - \v~u'\ 2H 


> \v 


-v' 


2H 


\u> 


-v\ - 


\v 


— 14/ 1 





(14) 
(15) 



(16) 
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by triangle inequality. Take u' — v = a and v' — u' = b. Then 

|« - v'\ 2H + \v! - v'\ 2H -\v- u'\ 2H = \a + b\ 2H - \a\ 2H - \b\ 2H (17) 

where < a < 1, < b < 1 and H G (1/2, 1). 
For < b < a < 1, one can obtain 

(a + b) 2H - a 2H - b 2H (18) 

= a 2H (l + -)-a 2H -b 2H 
a 

= a 2H (l + 2H- + 2g(2 f ~ 1} (1 + c) 2H - 2 (-f) - a 2H - b 2H 
a 2 a 

= b(2Ha 2H - 1 - b 2 "- 1 ) + H(2H - 1)(1 + c) 2H - 2 a 2H - 2 b 2 

by using Taylor expansion. Inequation ( fT8l , c e (0, b/a), 1 < 2i/ < 2 and < 6/a < 
1. Therefore, we see that 

(a + b) 2H - a 2H - b 2H = \v- v'\ 2H - \v! - v'\ 2H - \v - uf H > 0. (19) 

Then 

|u - u'\ 2H + \v~ v'f H -\u- v'\ 2H -\v- u\ 2H > 0. (20) 

For < a < b < 1, by symmetry in the calculations we get the same result. And for 
< a = b < 1, since 

(a + b) 2H - a 2H - b 2H = (2 2H - 2)a 2H > 0, (21) 

we obtain the same result. 

Hence, for case < U < v <«'<«'< 1, we see that; 

E{{BZ-B»)-{B%-B»)f > \u-v\ 2H + \u'-v'\ 2H 

> 2mm(\u-v\ 2H ,\u' -v'\ 2H ) 

2E(B*-B^) 2 , \u-v\ 2H < \u'-v'\ 2H 
2E(BH-B"f, otherwise 

Now repeating the same arguments in all 6 cases such as < u < v < u' < v' < 1, 
0<u<u'<v<v'<1,Q<u<u'<v'<v<1,0<u'<v'<u<v<1, 
< u' < u < v' < v < 1, < u' < u < v < v' < 1, one by one we see that some 
inequalities are satisfied again by the symmetry of calculations, that is 

E((B»-B»)-(B»-B»)) 2 

> 2 min(|u — v\ , \u — v \ , \u — v \ , \u — v\ , \v — v \ , \u — u \ ) 

> 2 min(|u — v\ 2 ,\u' — v'\ 2 ,\u —■ v'\ 2 ,\u' — v\ 2 ,\v — v'f ,\u — u'\ 2 ) 
= 2\s-t\ 2 

= E{V2B] - V2B}) 2 . (22) 
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Here, \s — t\ 2 = min(|u — v\ 2 , \u' — v'\ 2 , \u — v'\ 2 , \u' — v\ 2 , \v — v'\ 2 , \u — u'\ 2 ). 
Hence, by the Sudakov-Fernique inequality 

-^==E( sup V2Bl) = E(V2Sl) < E(M?) (23) 

V 7r 0<s<l 

is proved. □ 
Corollary 1. Consider fBm up to fixed time t with 1/2 < H < 1, 

^==E( sup V2Bl) < E(M t H ) < E{ sup y/2B l J 2 ) = 

V 77 0<s<t 0<s<t V 71 " 

and 

2t H 

P(M t H >x)< —=. 

Proof. By the self similarity property of fBm 

E(M t H ) = E( sup (B»-B»)) 

0<u<v<t 

= E( sup (t H (B* -B?))) 

0<u<v<l 

= t H E(M?). (24) 



And, by Markov's inequality, 



(25) 
□ 



2.2. Numerical Results 

Let Bf , t > be a fBm defined on the probability space (O, J 7 , P). We consider 
fBm with Hurst parameter 1/2 < H < 1, which is on the interval (0, t). We divide time 
t into n parts and obtain the increments {Bfg — B^_.^ s : i : 1, 2..., n; j = 0, 1, i} 
with 5 = t/n, where each one of them is centered Gaussian variable with different 
variances. After such discretization, we have the following inequality, 

max {Bg Bf t _ l)5 ) < M t H . (26) 

Covariance between increments The increments given above are centered Gaussian 
variables, and they are not independent, the covariance between them is 



E [( B i5 ~ B fi-j)s)( B kS ~ B fk-i)s)\ - 

+ 



E [ B iS B ks] - E [ B iS B fk-l)s\ 
E [ B (i-j)5 B ks\ + E \ B (i-j)S B (k-l)8\ 

±[\i5-(k-l)5\ 2H + \(i-j)5-k6\ 2H ] 
l -[-\i5-k5\ 2H -\{i-j)5-{k-l)5\ 2H ] 
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where i — 1, 2..., n; j = 0, l...,i;k = 1, 2..., n; / = 0, 1..., fc. 

In the literature, there are many studies on the expected value and the distribution 
of maximum of dependent Gaussian variables. The history of this issue starts with 
Tippett, lfl7l|. T eichroew. flql . Clark and Williams, [8] and Bose and Gupta, Owen 
and Steck, 1 1311 considered dependent but identically distributed Gaussian variables. 
Clark, J3l obtained first four moments for two dependent and differently distributed 
Gaussian variables, however these results were only for two variables. Ross, IU5I1 gives 
an upper bound for dependent and differently distributed Gaussian variable, however 
in our case we need a lower bound. Lai and Robinson, ifTl f ound another bound for 
dependent but identically distributed variables. Brodtkorb, [5] presented an approach 
for variables with singular covariance matrix, however this approach also is not proper 



for our study. There is no exact, theoretical solution for this problem. Ross, 1 14] 
presented numerical methods to obtain lower and upper bounds for the expected value 
of maximum of dependent and differently distributed (DDD) Gaussian variables using 
Theorem|2]given by Vitale, Qjj] . 

Theorem 2. Let Wi, Xi,Yi, i — 1,2..., n be DDD centered Gaussian variables. And 
for all i and j 

E((Wi - Wjf) < E((Xi - X,) 2 ) < E((Yi - Y 3 f) 
is satisfied. Then, for arbitrary constants rrii, i = l,2...,n, 

E[maxWi + mi] < £[maxl, + mi] < E^maxK; + mj. 



In our study, we implement Ross, ffl method in order to obtain lower bound for 
the expected value of maximum loss of fBm. Ross, 111411 numerically finds lower and 
upper bounds on the expected value of maximum of DDD Gaussian variables which 
can't be calculated explicity. Using covariances of DDD variables, perfectly dependent 
and differently distributed (PDDD) Gaussian variables and independent and differently 
distributed (IDD) Gaussian variables which satisfy Theorem [2] are found. Lastly, by 
Theorem [2] 

ElmaxiPDDD^] < E[max(DDD)] < E[max(PDDD 2 )] (27) 
E[max(IDDx)] < E[max(DDD)] < E[max(IDD 2 )] (28) 

are obtained. 
In our study, 

1. Let {Wi, i = 1,2..., n} be PDDD Gaussian variables and {B% - B%~ s , i = 
1, 2..., n, j — 0, 1..., i} are DDD Gaussian variables. Suppose 

E((Wi - Wi?) < E[((B» Bf*_ j)s ) (Si - Bf k _ l)S )f] 

holds. Then by Theorem|2]given above, 

E(max^) < E{wn{B5 - B^ s )) 
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One can construct many Wi's which satisfy the above inequality. Therefore, 
we try to find the maximum value of E(ma,xWi). The calculation of this is 

i 

given as E(maxWi) = E[E[maxWi\Z}} = J^° oo h(z)(j)(z)dz, where h(z) = 
E(maxWi\Z = z) and <fi(z) is the probability density function of standart nor- 

i 

mal variable. 

2. Let {Yi, i — l,2...,n} be IDD Gaussian variables and {B?$ — B^_^ g , i = 
l,2...,n,j = 0,1..., i} are DDD Gaussian variables. Suppose following in- 
equality holds, 

Em - Yj) 2 ) < E[((B* - B%_ j)s ) (B* Bf k _ l)5 )f] 
Again, by Theorem[2] we have 

EipaxYi) < E(mzx(B? s - Bf^ j). 

Hence, same algorithm as the algorithm given in 1. is applied. The only differ- 
ence is, in calculating E{maxYi) that is, 

i 

^(maxr,) = / (1 - TT Pr(Y t < y) - TT Pr(y, < -y))dy 
J o i=l i=l 

Following this algorithm, we have conducted a numerical study using MATLAB. Table 
Q] shows numerical lower bounds we have found on the expected value of maximum 
loss of fBm up to time 1, with Hurst parameter H € (1/2, 1). 







n = 5 


n= 10 


n = 20 


n = 30 


n = 40 


n = 50 


Hi 
o 

II 


IDD Lower Bound 


0.6350 


0.6551 


0.5753 


0.5171 


0.4753 


0.4434 


tq 


PDDD Lower Bound 


0.2847 


0.1621 


0.1294 


0.1056 


0.0903 


0.0872 


CO 

o 

II 


IDD Lower Bound 


0.5230 


0.4855 


0.3954 


0.3407 


0.3040 


0.2772 


m 


PDDD Lower Bound 


0.2429 


0.1250 


0.0956 


0.0579 


0.0624 


0.0535 


O 

II 


IDD Lower Bound 


0.4309 


0.3550 


0.2645 


0.2179 


0.1886 


0.1680 


tq 


PDDD Lower Bound 


0.2037 


0.1055 


0.0695 


0.0491 


0.0425 


0.0305 


00 

o 

II 


IDD Lower Bound 


0.3578 


0.2688 


0.1830 


0.1432 


0.1196 


0.1038 




PDDD Lower Bound 


0.1682 


0.0859 


0.0478 


0.0357 


0.0274 


0.0192 


°i 
o 

II 


IDD Lower Bound 


0.2760 


0.2010 


0.1324 


0.1005 


0.0819 


0.0696 


tq 


PDDD Lower Bound 


0.1420 


0.0674 


0.0348 


0.0201 


0.0200 


0.0127 



Table 1 : Numerical lower bounds on the expected value of maximum loss of fBm up to time 1 
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2.3. Simulation Results 

In our study, we generate fBm with l/2<iJ<lupto time 10000 with 10000 
simulations using Cholesky method. This method is an exact simulation method for 
fBm, iflQIl . From this, we collect 10000 values of maximum loss of fBm up to time 
10000. We find their average, which gives us the expected value of maximum loss 
of fBm up to time 10000. Afterwards, using the self similarity property of fBm we 
calculate the expected value of maximum loss of fBm up to time 1, as follows 

E (M* ) = (29) 
We repeat this for different values of Hurst parameter. The results are given in Table|2] 



Hurst Parameter Cholesky Simulation results 



0.5 1.239 

0.6 1.00721 

0.7 0.82509 

0.8 0.69865 

0.9 0.61016 



Table 2: Simulation results for the expected value of maximum loss of fBm up to time 1 



3. Conclusion 

Any information on maximum possible loss as a measure of risk would be impor- 
tant for investors in order to hedge it or manage it. Here, we provide an upper bound 
and a lower bound on the expected value of maximum loss of fBm up to time 1 and up 
to fixed time, t. By using Markov's inequality, we find an upper bound on the distribu- 
tion of maximum loss. Later, we obtain numerical lower bounds on the expected value 
of maximum loss of fBm up to time 1 with 1/2 < H < 1. Table [3] shows all these 
results on the expected value of the maximum loss of fBm up to time 1 . 



Hurst Cholesky Simulation Caglar and Vardar Caglar and Vardar Vardar and Cakar Vardar and Cakar 

Parameter Results Lower Bound Upper Bound Lower Bound Upper Bound 

(1.5 1.239 0.39904 1.59617 0.56418 1.12866 

0.6 1.00721 0.39904 1.59617 0.56418 1.12866 

0.7 0.82509 0.39904 1.59617 0.56418 1.12866 

0.8 0.69865 0.39904 1.59617 0.56418 1.12866 

0.9 0.61016 0.39904 1.59617 0.56418 1.12866 

Hurst IDD Lower Bound PDDD Lower Bound 

Parameter n-5 n - 10 » = 20 n = 30 n = 40 n = 50 n =5 n = 10 n = 20 n =30 n = 40 n = 50 



0.5 0.6350 0.6551 0.5753 0.5171 0.4753 0.4434 0.2847 0.1621 0.1294 0.1056 0.0903 0.0872 

0.6 0.5230 0.4855 0.3954 0.3407 0.3040 0.2772 0.2429 0.1250 0.0956 0.0579 0.0624 0.0535 

0.7 0.4309 0.3550 0.2645 0.2179 0.1886 0.1680 0.2037 0.1055 0.0695 0.0491 0.0425 0.0305 

0.8 0.3578 0.2688 0.1830 0.1432 0.1196 0.1038 0.1682 0.0859 0.0478 0.0357 0.0274 0.0192 

0.9 0.2760 0.2010 0.1324 0.1005 0.0819 0.0696 0.1420 0.0074 0.0348 0.0201 0.0200 0.0127 

Table 3: Results for the expected value of the maximum loss of fBm up to time one 
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Here, Caglar and Vardar upper and lower bounds are the theoretical bounds given 
in ffj. Comparing the simulation results with the values given in Table [3] we see that 
Vardar and Cakar theoretical upper bound and lower bound which were obtained in 
TheoremQ] for 1/2 < H < 1, provides the closest upper and lower bounds among all 
the bounds. 
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